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The so called f{X) hybrid metric-Palatini gravity presents a unique viable generalisation of the 
f{R) theories within the metric-affine formalism. Here the cosmology of the f{X) theories is stud¬ 
ied using the dynamical system approach. The method consists of formulating the propagation 
equation in terms of suitable (expansion-normalised) variables as an autonomous system. The fixed 
points of the system then represent exact cosmological solutions described by power-law or de Sitter 
expansion. The formalism is applied to two classes of f{X) models, revealing both standard cos¬ 
mological fixed points and new accelerating solutions that can be attractors in the phase space. In 
addition, the fixed point with vanishing expansion rate are considered with special care in order to 
characterise the stability of Einstein static spaces and bouncing solutions. 

PACS numbers: 


I. INTRODUCTION 

The so called hybrid f{X) class of modified theories of gravity consists of the superposition of the metric Einstein- 
Hilbert Lagrangian R with an f{TZ) term constructed a la Palatini [l|. This interpolation between the non-equivalent 
metric and Palatini approaches has revealed interesting possibilities. For instance, in the formulation of the dynami¬ 
cally equivalent scalar-tensor representation, it was shown that the theories contain a “light”, long-range scalar field, 
which is able to modify the cosmological and galactic dynamics while leaving the Solar System unaffected [li- This 
feature has led to investigations of the cosmology in f{X) theories P-Q, but a systematic mapping of its phase space 
is still missing. It is this gap we aim to hll in the paper at hand. 

From a theoretical point of view, the f{X) class of theories enjoys a similar special status amongst the more general 
hybrid metric-Palatini theories as the f{R) theories within the narrower framework of purely metric gravity. To wit, 
when one excludes theories which present ghost-like, superluminally propagating and otherwise pathological degrees 
of freedom, there is evidence^ that the f{X) class of theories is singled out as the only viable form of an action 
one can construct using the metric, its derivatives (and thus the metric Levi-Civita connection) and an independent 
’’Palatini connection” [T^ . The underlying reason is [HI that in the special case of f{X) actions the higher derivatives 
in the gravity sector, that otherwise would result in an Ostrogradskian instability, can be separated into a scalar 
mode. Indeed the f{X) gravity, where X = R -\- T, represents a generic case within the one-parameter family of the 
Algebraic Scalar-Tensor theories, at one end of which lies the pure Palatini f{TZ) (wherein the field is a function of the 
stress energy trace T) and at the other end the pure metric f{R) (where the field is a function of the metric curvature 
J^) |l2l|. Furthermore, the propagating degrees of freedom have proven to be healthy also on curved backgrounds 
[3|, l7| , and concerning the Cauchy problem, it was shown that in this class of theories the initial value problem can 
always be well-formulated and well-posed depending on the adopted matter sources [l^ . These considerations provide 
compelling motivation for further exploration of these particular theories. 

An immediate application is to consider gravitational alternatives to dark energy. As highlighted above, the 
hybrid theories are promising in this respect as they can naturally avoid the local gravity constraints but modify 
the cosmological dynamics at large scales. This is simply because as a scalar-tensor theory, the hybrid f{X) gravity 


* Electronic address: sante.carloni@tecnico.ulisboa.pt 
"I" Electronic address: tomi.koivisto@nordita.org 
^ Electronic address: fslobo@fc.ul.pt 

^ In particular, we can exclude the more general cases of f{R, 'IZ) studied in and various second, fourth and sixth order derivative gr avity 
theories resulting from considering more general curvature invariants constructed from the metric and the two available connections [ig| . 



2 


is characterised by an evolving Brans-Dicke coupling^, which allows to introduce potentially large deviations from 
General relativity in the past (and future) as long as the coupling at the present epoch is strong enough to hide the 
field from the local gravity experiments. In contrast, in the metric f{R) models the Brans-Dicke coupling is a finite 
constant and one needs to invoke some of the various ’’screening mechanisms” 0 (workings of which remain to be 
studied in the hybrid theories) in order to reconcile the Solar System experiments with cosmology. Such dark energy 
dynamics that drive the theory towards its general relativistic limits today have indeed been realised in several specific 
models A preliminary phase space analysis confirmed the existence of such de Sitter attractor solutions, and 

also other analytic solutions were presented in Ref. Q as well as in Ref. Q, there using a Nother symmetry technique. 
A designer approach was employed in Ref.Q to reconstruct precisely the standard ACDM expansion history by a 
nontrivial f{X) model, and finally, two families of models were constrained by confronting their predictions with a 
combination of cosmic microwave background, supernovae la and baryonic accoustic oscillations background data Q . 
At this point, it can be useful to have a global view of the possible cosmological solutions from the perspective of 
dynamical system analysis. 

Cosmological perturbations have been also analysed in these models up to the linear order , and the results 

imply that the formation of large-scale in the aforementioned dark energy cosmologies is viable though exhibits subtle 
features that might be detectable in future experiments. The full perturbations equations were derived and at their 
Newtonian limit, describing the observable scales of the matter power spectrum, the growth of matter overdensities 
was shown to be modified by a time-dependent effective fifth force that will modify the red-shift evolution of the 
growth rate of perturbations 0,i. In addition, numerical studies of the perturbations imply that the difference of 
the gravitational potentials can exhibit oscillations at higher redshifts even when the background expansion and the 
full lensing potential are indistinguishable from the standard ACDM predictions Q. Such features could potentially 
be observed in cross-correlations of the matter and lensing power spectra, but detailed calculations of the cosmic 
microwave background anisotropies and other perturbation observables remain to be carried out. 

At an effective level, the f{X) modifications involve both (the trace of) the matter stress energy and (the Ricci 
scalar of) the curvature, and from this point of view it appears appealing to speculate on the possible relevance 
of these theories to both the problems of dark energy and dark matter, in a unified theoretical framework and 
without distinguishing a priori matter and geometric sources [^. Various aspects of dark matter phenomenology from 
astronomical to galactic and extragalatic scales have been considered in the presence of f{X) corrections to gravity 
with some promising results ii Il5| . The generalised virial theoreom can acquire, in addition to the contribution 
from the baryonic masses, effective contributions of geometrical origin to the total gravitational potential energy, 
which may account for the well-known virial theorem mass discrepancy in clusters of galaxies (H. In the context of 
galactic rotation curves, the scalar-field modified relations between the various physical quantities such as tangential 
velocities of test particles around galaxies, Doppler frequency shifts and stellar dispersion velocities were derived @ . 
More recently, observational data of stellar motion near the Galactic centre was cornmred with simulations of the 
hybrid gravity theory, which turned out particularly suitable to model star dynamics [l5j. Yet, to promote the f{X) 
theory into a convincing alternative to particle dark matter, one should also produce the cosmological successes of 
the ACDM model without the CDM component. 

In this paper we study the cosmological solutions of the theory in terms of dynamical system analysis. Cosmological 
applications of the phase space techniques have a long history, dating back to at least the early studies of chaos in the 
framework of cosmology by Belinsky (l^ and Misner 0, and were developed in particular by Ellis and Wainwright 
(see [3 references therein). The nowadays standard framework of constructing the phase space is based on 
defining expansion-normalised variables representing the various (effective) contributions to the total energy density. 
Besides mathematical efficiency this approach allows an easy physical interpretation of the phase space features as 
its variables assume the structure of cosmological parameters. The method has been used with great success in the 
analysis of anisotropic cosmologies [l^ . and has been applied e.g. to models with nonstandard equation of state 0, 
and more recently exploited to analyse a number of modified gravity models^ In the context of f{X) gravity, 
the stability of the Einstein universe, a particular static fixed point in the phase space of the theory, was analysed 
in great detail using the equivalent scalar-tensor description |22l |. and matter-dominated and de Sitter fixed points 
have been checked also Here we will set-up the phase space directly in terms of the curvature quantities, with 
the immediate advantage of a more direct interpretation of the relation between the cosmic phenomenology and the 
form of the fundamental function f{TZ), and moreover allowing to take full advantange of the conformal relation 


^ To be precise, we have uibd = 3<?i/(2<^ — 2Q,a), where = 1 for our case. This normalisation is chosen to produce the observed value 
of the Newton’s constant today. The pure Palatini case is > 0, and the purely metric limit is recovered at fix —f oo. 

® To highlight the power of the method, we could mention that subjected to its scrutiny, the rather involved string-inspired set-up in 
an anisotropic spacetime filled by a massive vector field and a scalar field, and featuring their both conformal and disformal couplings, 
yielded to complete and detailed uncovering of its phase space l2lll . 
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between the metric and Palatini curvatures and resulting in a more efficient and generally applicable formulation of 
the dynamical system. 

This paper is outlined in the following manner: In Section m we present the basic formalism for the hybrid f{X) 
gravity, namely, the action and the field equations and adapt them to cosmology. In Section uni we set up the 
dynamical system, which essentially consists of defining suitable variables and deriving their propagation equations. 
We then analyse the phase space structure in several families of particular models. In Section El we reformulate 
the system in order to analyse the fixed points with H = 0, that reside in the asymptotic regime of the phase space 
considered in the previous section. In Section |Vl we discuss our results and conclude. 

As the reader may have noticed, we have set 8ttG = 1. Natural units {h = c = ks = 8 ttG = 1) will be used 
throughout this paper, and Greek indices run from 0 to 3. We use the (—,+,+,+) signature and the Riemann tensor 
is defined by 

uP — a r 0 — a r ^j-r ^ — P ;5^P ^ ('ll 

uoc G/j/x \ fid ^ ya. ^ yd ^ fia. 5 V-*'/ 

where the are the coefficients of a linear dynamical connection. The Ricci tensor is obtained by contracting the 
first and the third indices via the metric i.e., = R°‘^av- 


II. BASIC EQUATIONS 

The four dimensional action^ of hybrid metric-Palatini gravity is given by [ij 

5 = i y d^xv^ [R + /(7^) + V')] , 


( 2 ) 


where Lm = is the matter Lagrangian, coupling tires matter fields ip minimally to the metric R is 

the Einstein-Hilbert term, TZ = g^'^TZ^n, is the ’’Palatini curvature”, given by the contraction of the tensor that is 
generated by an a priori independent connection P“^, as 


_ "pa _ pa I pa pA _ pa pA 

''^fiy — -L ^ fiOL,y ' ^ aA"^ fiy ^ fi\^ cty ' 


Varying the action ([2]) with respect to the metric, one obtains the following gravitational field equation 


Gfiv + F{TZ)TZfj_i, — —J(TZ)g^v — T/xi/; 


with the notation F{TZ) = f'ijV) and the usual definition of the matter stress-energy tensor 

rp _ 2 5{y/—gLm) 


(3) 


(4) 


(5) 


The field equations for the connection imply that the independent connection is compatible with the conformal metric 
g^i, = F{TV)g^i, and the relation 




2F'^{ny 


F{n) 


2F{ny 


( 6 ) 


connects the Palatini and the metric Ricci tensors. 

Here we will focus on the cosmology of the hybrid models. We consider the homogenous and isotropic cosmologies 
represented by the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric 


ds^ = —dt^ -I- a(t)^ 


dr^ 


1 — fcr^ 


(d6|2 sin^ 


(7) 


^ In this setting / contains a coupling constant, but since it will be integrated in one of the variables, its presence will be irrelevant for 
the following discussion. 
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The Hubble rate H of the scale factor a is defined as H = a/a, where an overdot denotes a derivative with respect to 
the time coordinate t. The Friedmann equations for this Hubble rate, derived from the above general equations, can 
be now written as Q 

Peff = - 3fc, (8) 

Pee = -2H - 3H^ + 3k , (9) 


where peff and Peff are effective energy density and pressure given by 


Pee = p-\\fiX)-^ 


F"{X) - 


{F'{x)y 


Fix) J 


- 


^F'(X) (x + XH^ - 3F{X) (h + } 


Pee = P + I i/(X) - iF"(X)X2 - iF'(X) (x + 5Xh) - F{X) (h + 3H^) } , 


( 10 ) 

( 11 ) 


respectively, where p and p are the energy density and pressure of a perfect fluid, and X is defined as 

X = T + i? = F(7^)7^ - 2/(7^), (12) 

where we have used the trace of Eq.(|3]) in the last equality. Note that in the pure Einstein gravity one has X = 0, and 
the curly bracket terms simply vanish: this formulation of the cosmological equations has the advantage of making 
explicit the deviation of this class of theories from General Relativity and placing in evidence the scalar degree of 
freedom introduced by the Palatini term of the action ([2]). 

As the next step, we take advantage of the conformal relation between the two geometries ([5]) that naturally appear 
in the theory. Indeed, the equations (|5]) and ([51) can be simplified with a suitable choice of variables that exploits in 
each occasion the optimal quantity associated with either or g^y. Defining the auxiliary function A = yjF{TZ) a(t), 
one can write 


n 

n 


a F _ A 
a ^ “ A’ 




(13) 

(14) 


where the symbol f indicates the derivative with respect to r, defined by dr = ^FijVjdt. Using these quantities 
defined in Eas. (fT^ - (IT^ to rewrite the field equations (|4|), they assume the following form: 

f-) +^il + F) + F'H^-\iFn-= (15) 

\a J D o 

- — H— ^ ^ (^ ”*" ~ ® ’ (^®( 

a \ a / DO 

fi + 3^(1 + w)p = 0 , (17) 

where we have used the equation of state p = wp, and for notational simplicity, we have dropped the dependence 
of / and F of TZ. This formulation of Eqs. m and ([2]) can be used to deduce some basic features of the hybrid 
metric-Palatini theories. For example, Eq. (1161) shows clearly that the interaction between the Palatini and metric 
part of the theory can generate an increase of the expansion rate for a flat and spatially closed universe {F must 
be definite positive for the transformation g^y = FiTVjg^y to make sense), whereas the nonminimal coupling of the 
Palatini curvature can have the same effect only if the function / is negative. This means that a sufficient condition 
for the model to generate cosmic acceleration is / < 0 and k = 0,1 (at least when the geometrical contribution 
dominates the matter sources). Similar considerations can be used in other frameworks such as in the context of (a 
homogenous and isotropic) gravitational collapse. 

In the following we will use the above equations to implement the dynamical system analysis. Differently from 
previous similar attempts, we will not rely on the definition of an auxiliary scalar field to perform this task, but we 
will only use Eqs. (fT5]l - (frh . 









5 


III. GENERAL DYNAMICAL SYSTEMS FORMALISM 

In this section, we will construct the dynamical system for the hybrid models, which consists of essentially two 
steps: (i) the definition of suitable variables, and (ii) the derivation of the propagation equations for these variables. 
Basically this is rewriting of the system of differential equations as a first order system of the chosen different variables. 
The fixed points will thus correspond to particular exact solutions of our original Friedmann equations (|5]) and ([5|). 


A. General formalism 


Let us define the variables 


A = —, y = ^, z=^, n = K=^^. 


( 18 ) 


The interpretation of these variables is straightforward, as each represents the corresponding quantity in units of the 
expansion rate. We also introduce a dimensionless time variable: the logarithmic time A/” = Ina. Note that with this 
time variable we are implicitly imposing that the orbit of the phase space will represent only monotonically expanding 
or contracting cosmologies. We will focus on the possibility of solutions crossing iL = 0 in the following section. 

In terms of the variables (ITSl) , the cosmological equations can now be written as the following dynamical system: 


A At = [2Z - X - 2T {K + A^) + (1 + 3w)il] - K+ Y , 

Y,j^ = Y [2(1 + Z)-2T{K + A2) + (1 + 3w)n + 2(A - 1)Q] , 

Z At = Z [2(1 + Z)-2T{K + X^) + (1 + 3w)n] + 2Y{X - 1)AQ, (1^) 

0,At = ri[2Z- 2T {K + A^) + (1 + Zw){n - 1)] , 

A At = A [2Z - 2F (A + A^) + (1 + 3w)ff] . 


with the constraint 


A[1 + A] + A^A -YA-ff + Z + l = 0. 


( 20 ) 


In the above equations, the quantities A = F{TV) and Q(A) = F/{TZF') are only functions of TZ, and in 
the system they have to be expressed in terms of the variables specified in Eq (fT^ . Given the form of 
done by noticing that Z/Y is function of TZ only, i.e., 

^ /( 7 ^) 

Y TZ ■ 


order to close 
/, this can be 


( 21 ) 


Inverting this relation for Y ^ 0, one obtains TZ = A(Z/Y) and the dynamical system can be closed. Therefore, as 
far as this inversion is possible, the phase space spanned by our dynamical variables (|18D is fully described by the 
closed system of propagation equations (1191) . The price to pay for this is the same that one finds in the application 
of a similar method to the metric f{R) gravity |^. 1^. namely that the phase space can have singular surfaces (such 
as Y = 0), which may have to be treated with great care. 

In the following we will use the constraint equation (1201) to eliminate the variable If. The independent equations of 
the system are 


A aa = ^X {A [(3w - 1) (A + A2) - (3w + 1)Y] + (1 + 3w)(K + 1) + 3(w + 1)Z - 2A} - A + Y, 

Y;a^ = Y {A [{3w - 1) (A + A^) - {3w + 1)Y] + (1 + 3w)K + 3{w + 1)(Z + 1) + 2Q(A - 1) + 2} , ^^2) 

Z^^f = Z {A [{3w - 1) (A + A2) - (3w + 1)Y] + (1 + 3w)K + 3{w + 1)(Z + 1) + 2} + 2AQY(A - 1), 

A aa = a {a [{3w - 1) (a + A^) - {3w + 1)Y] + (1 + 3w){K + 1) + 3{w + 1)Z} . 

Any other choice of independent dynamical variables could be made in principle, but the one above is the only one 
that does not introduce divergences in the dynamical equations. 

The system above can then be analyzed with the standard tools of dynamical system analysis and in particular 
with the Hartmann-Grobmaim theorem. Once the fixed points have been found, together with their stability, the 
solutions associated to the fixed points correspond to 

a = ao{t - , (23) 

p = (24) 
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where the constant 7 is obtained by expressing the modified Raychaudhuri equation (1161) in terms of the dynamical 
variables (the asterisk denoting their respective value at the fixed point) 

7 = + XlT^ - ^(1 + (25) 

If 7 is equal to zero the formulae above do not hold and we obtain an exponential solution, which can be fully specified 
as well by substitution into the cosmological equations. 


B. Specific case I: / = yT?." 

As an example let us look at the specific case / = We have, for R 7 ^ 0, 

f = 


( 26 ) 


so that 




(27) 


It is important to point out here that in the case x = 0 fhe (I26p cannot be inverted. This result implies that we will 
always consider x 7 ^ 0 in our analysis®. The system (1^^ reduces to 


[(3w- 1) (K + X^) - (3u; + l)y] + {1 + 3w){K + 1) + 3{w + 1)Z - 2x'^ - K + Y, 


YM = y 


(1 + 3w)K + 3 ( 1 ^ + ^){Z + 1) + 2--—h 2 

n — 1 


+ Z [(3t/; -1){K + X‘^) - {3w + 1)Y] 


Z j\f — ^ [(Suj — 1) {K + — {3w +(1 + 3w)K H- 3{w H- 1) (Z H- 1) + 2 H - — [X — ^ ’ 

= ix [( 3 iy - 1) (a: + X^) - {3w + l)y] + (1 + 3w){K + l) + 3{w + l)z| . 


(28) 


It is important to stress at this point that the system above is equivalent to only if T 7 ^ 0, so fixed points with 
y = 0 are excluded. This fact remains true even if no denominator containing Y is present. The reason is that the 
(l?T]) only holds for T 7 ^ 0 and therefore for F = 0 one cannot close the system. Note also that the (1^51) now includes 
two invariant submanifolds, Z = 0 and K = 0. This result implies that only points with Z = 0, AT = 0 (and fl = 0) 
can be global attractors for the phase space. 

The fixed points are determined setting to zero the A/'-derivatives of the dynamical system variables. The finite 
phase space contains three fixed points. Their associated solutions and stability can be found in Table HI We can 
distinguish two classes of fixed points depending on the value of the variable Z. Points with Z = 0 {B,C) essentially 
represent states in which f{TZ) ~ 0 and therefore in which this model approximates General Relativity. 

It is also worth noting that the points B and C, can be characterized by a negative value of X i.e. has the peculiar 
characteristic to have negative "Palatini Hubble rate" while the physical universe is expanding or vice versa. 

All the solutions associated to the fixed points are the typical GR Friedmann solutions already found in ii- In 
addition, the coordinates of the associated fixed points reveal that GR-like solutions can only be realized Z = 0, as 
expected. Point A is associated to an exponential expansion whose time parameter A can be obtained from the trace 
of the gravitational field equations (fT^ . In the case of / = TZ^ this equation reduces to 


(n - 2)7^ = T + i?. 


(29) 


Then, in the case of an exponential solution in vacuum, R is constant and 


7 ^ = 


R 

n — 2 


(30) 


^ This conclusion is in general true for every type of if the coupling constant in / is zero our entire construction cannot be applied. 
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Table I: The fixed points and the solutions of the hybrid model with / = yTT.". Here the value of lambda is given by 
the (l?T]) . In addition, A stays for attractor, for attractive focus, R for repeller and S for saddle. 


Point (A, y, y, U, K) 

Scale Factor 

Energy Density 

Stability 

A {l,2,;^,0,0} 

a = flo exp [X{t — to)] 

p = 0 

Fa 

B {2-n,(n-l)(n-3),0,0,-1} 

O 

1 

O 

p = po (t - to) 

Fa 1 < n < 2 

0 < w < 1 

S otherwise 

{|[-3n(M) + 1) + 3u; + 5], 

C j[3n{w + 1) — 4][3n(ui + 1) — 3w — 
0,1,0} 

5], a = ao{t — to) 3(i+“) 

p = Po(t — to)~^ 

S 


This immediately implies that H =% and substituting into the Friedmann equation one obtains A = 0 and 

1 


A = 


2^3 


(n - 2)2("-l) [(4 - n)n - 2]2~2n^2-2n . 


(31) 


For a given n this expression will give multiple solutions. This is a common occurrence in higher order gravity and it 
is also present in the hybrid metric-Palatini case. However, the problem arises of how to associate these values of A 
to the solutions in the fixed points. The case A < 0 can be excluded if we consider only expanding cosmologies, and 
the case A = 0 because it implies H = 0 that cannot be considered finite phase space. Further selection criteria can 
be obtained by checking the values of the dynamical variables at the fixed points. In the case of A, one has T > 0 
which implies TZ > 0. This result allows to exclude all the values of A associated to 7?. < 0. Therefore one is left with 
the (unique) real positive value of Eq. (|31ll which requires TZ > 0. Unfortunately, no general expression for this value 
of A can be given explicitly, but its calculation is straightforward for a specific n. 

Looking at the stability of the fixed points, there are only two possible attractors: A and B. The last point behave 
as a sink only for specific values of n. Since, however, these points do not have Z = 0, K = 0 they cannot be global 
attractors. In fact, in general no fixed point can be a global attractor for the cosmology of this model. The points 
with Z = 0 are generally unstable for every value of the parameter n. The exception is given by the point B for 
1 < n < 2. This result implies that GR-like Friedmann cosmologies are in general unstable in this model and only a 
special set of initial conditions might lead to the (non global) attractor B. Such feature can be interesting to model 
the onset of a dark era in this context. 


C. Specific case II: / = 

As a further example, let us consider the case of the function / = aTZ/{l + PIZ), with a and (3 generic constants. 
This model is inspired by the proposal of Hu and Sawicki [2^ in the context of the metric f{R) gravity and we will 
refer to it as “hybrid HS” model. As in that case the function / is such that in expanding cosmologies for very high 
and very small values of the Hilbert-Einstein term R, T is also very large or very small and the Palatini curvature 
TZ is constant. In particular with our choice of the values of the parameter, TZ becomes zero (a positive) or — ^ (a 
negative) for R —0 and becomes — ^ if R —>■ c». Since we also require that F > 0 the constant a must be chosen 
positive. Thus effectively the function / represents an effective cosmological constant term which appears in the late 
universe, but does not affect the early one. Note that the sign of this effective cosmological constant depends on the 
sign of /?, so that if /3 < 0 the effective cosmological constant in Eq. O is positive otherwise, if /3 > 0, it is negative. 
We will choose, in the following, the first option. 

It is also worth noticing that for this type of function /, in the limit TZ —>■ (0, oo), F becomes constant and Ti = H 
so that the equations effectively reduce to GR plus a cosmological constant. This is a characteristic property of the 
Hu-Sawicki model in the context of metric /(R)-gravity and it is present also in the hybrid metric Palatini case. 









For the hybrid HS model we have, for V 0, 


so that 


Z a 

r ~ 1 + /37Z 


Z^ 


Q 


aY 

2{Z -aY)’ 


(32) 


(33) 


Note that in the case /3 = 0, f{TZ) becomes linear. This case is degenerate in terms of the dynamical system variables 
{Y = Z) and has to be treated separately®. Using (1551) . the system (1551) becomes 


Xx = - 1) (^ + - (3w + 1)U] + (1 + 3w)iK + 1) + 3{w + 1)Z - 2x| - K+ Y 


Z^ 

Y^ = — [(3t„ -1){K + X^) - {3w + l)r] + Y 


Zx — z 


(1 + ?,w)K + i{w + \){Z + 1) + y + 2 

Z — ai 

Z{X-1)- 


{ [(3w^ — ^) {K + X"^) — {3w + 1)F] + (1 + Sw)K + 3{w + 1){Z + 1) + J | 


(34) 


Kx = [(3«^ - 1) (^ + - (3w + 1)U] + (1 + 3w){K + 1) + 3{w + l)z| . 


The above system only admits the invariant submanifolds Z = 0, A" = 0 and presents a singular subspace for Z = aY. 
The analysis of the phase space will be performed in the Z > aY and Z < aY part of the space with the regular 
techniques. We do not address here the behaviour of the orbits close to Z = aY. 

The system (1551) admits four fixed points as indicated in Table El Of these points, two {A, B) correspond to de 
Sitter solutions and the other two have Z = 0 and represent GR-like solutions. 

Again, the time constant for the de Sitter solutions can be calculated using the trace equation, rewritten in Eq. 
dm)- However, for this case the solution for the Palatini curvature is given by 


a + 24/3A^ =F 

4/3 (a + 62A2) 


(35) 


which is not unique. Therefore, differently from the previous case and the f{R) case, here the de Sitter solutions are 
set apart by the value of the Palatini curvature, other than R. As in the previous example the specific value of A 
associated to a fixed point can be deduced via considerations on the type of evolution represented by the phase space 
and the coordinates of the fixed points. In particular, A = 0 and A < 0 as well as the values of A for which TZ diverges 
have to be discarded. This leaves only two values of A, given by 


Ai 


1 ^ I 1 + a + y/a{a + 1) 

2Y 59 


A2 


1 ^ Ia + 1 — ^/a{a + 

2Y 59 


(36) 


respectively. Assuming /3 < 0, we have that, for 0 < a < 1/3, Ai is associated to 7^ > 0 and A 2 to 7^ < 0. Instead, 
for a > 1/3, Ai is associated to 7?. < 0 and A 2 to 7?. > 0. As a consequence, the nature of the fixed points depends on 
the value of a other than their coordinates. Specifically, points A and B that have F > 0 are characterised by Ai for 
0 < a < 1/3 and A 2 for a > 1/3. 

The calculation for the stability for these fixed points is straightforward. The points A and B are always stable. 
The phase space present no other possible attractor. Since global attractors can only be flat GR solutions (Z = 0, 
AT = 0), in spite of the generation of an effective cosmological constant discussed in the beginning, our analysis shows 
set of initial conditions for which a de Sitter expansion is not realised. In other words, the phenomenon of approaching 
to a future dominance of an effective cosmological term is not obvious as our initial reasoning seems to indicate. An 
escape from this conclusion could be the presence of some asymptotic {H = 0) attractor. We will see in the following 
that this is not the case. 

Notice also that in this model there is no GR-Friedmann past attractor. However considering also the stability of 
point V which is related to the classical Friedmann solution we can conclude that GR-like Friedmann cosmologies are 
always unstable. 


® In fact upon closer inspection, it appears clear that if = const, then A oc a and 'H = H and the dynamics is governed only by the 
variable ft. The phase space is one dimensional and presents only a single, unstable, fixed point corresponding to a t2/3(l+ui) solution. 





















Table II: The fixed points and the solutions of the case with / = The value of the parameter A is 

determined by Eq. (1551) . Here A stands for an attractor, R for a repeller, S for a saddle, for an attractive focus 
and F/j for a repulsive focus. The stability of point J is given only for the case w = 0. 


Point 

{X,Y,Z,n,K) 

Scale Factor 

Energy Density 

Stability 

A 

[l,2,2a-2^a(Q + l),0,0] 

a = ao exp [\{t — to)] 

p = 0 

Fa 0 < a < 49/32 
A otherwise 

^ 1 

l,2,2o + 27 a(a + l),0,0] 

a = ao exp [A(t — to)] 

fj, = 0 

Fa 

C 

[3,8,0, 0,-1] 

a = ao {t — to) 

fj, = flo{t — to)”^*'^'*''"^ 

S 

V [4 + 3w,i(3w + 4)(3«; + 7),0,1,0] 

a = no (t — to) 3(i+”) 

p = po(t — to)~^ 

S 


IV. FIXED POINTS WITH H = 0 


The dynamical system formulation above is not the only one that can be constructed for the cosmological equations 
(HHl-IIll). Another possibility is to use the following variables 


X=^ Y=^ Z=^ 


H = 




K = 


kF 


(37) 


where T-L = ^JF{^VjA^/A = A/A, together with the time variable A1 defined in such a way that, for a given dynamical 
system variable P, Pm = The previous set of variables differs from this case, as here one is able to see fixed 

points characterized by = 0 which are in the asymptotic regime of the previous phase space. Using the (1371) we are 
able to characterize the stability of Einstein static universes in the hybrid models and to gain information on bounce 
{H = 0, H > 0) phenomenology or other changes in the sign of FI typical of higher order theories of gravity (^ . 

Before proceeding with the analysis, it is important to stress that the new time coordinate M differs in sign from 
Ai whenever A < 0. This fact implies that in treating the stability, fixed points for which A < 0 will have different 
stability (e.g. attractors could correspond to repellers for given values of the parameters) with respect to the ones 
obtained with the variables of the previous section. As we will see in the example we considered the stability of the 
fixed points does not change, but this might be the case in other modelts. It is clear however that, since H is related 
to the physical time the “correct” stability is the one obtained in the previous section. 


A. Dynamical system equations 


The dynamical system equations in this case read 


X^M = {K+ 1)F + X{K - r + 1) - i(l + 2,w)n - Z , 

(38) 

Ym = 2U[i? + A - F - (A - 1)Q + 1], 

(39) 

Z^M = ‘2Z[{K + A- F + 1)-(A- l)YFQ , 

(40) 

- (3w + 1)A - 2r + 2], 

(41) 

K,m = 2A(1 + K-Y), 

(42) 

with the constraint 


(1 + A - Y)F + A + A2-H + Z = 0, 

(43) 

and the evolution equation for % 
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which is decoupled from the other dynamical equations. 

Analogously, in the following, we will eliminate the variable O for consistency. In this way the independent equations 
of the system are the following 


X.M=X 


K-Y 




Ym = ‘^Y[K + X-Y-{X-1)Q+1], 

Z^M = ^[{K + X-Y + 1)Z-{X- l)YFQ], 
K^m = 2K{1 + K-Y). 


(45) 

(46) 

(47) 

(48) 


The solutions associated to the fixed points in this case are given by equation is Eq. ((4^ rather than the Raychaudhuri 
equation (which is used only when 'H{t) has been determined in the fixed point). We have 


H = - 


X^ + Z* - (1 + + -(1 + 3w)n^ 


n. 


(49) 


where the asterisk means, as usual, that the variable is evaluated in the fixed point. From this result one can derive 
a and fi. 

In the following we will repeat the analysis of the / = yTT." and the hybrid HS model using these variables. 


B. The case / = yTT." in the barred variables 


In this case the general system reduces to 


X,M=X 


K-Y+ l-^{l + 3w)X 


1 nZ 

--[K{1 + 3w) + 3(1 + w)Z - —[(I - 3w){K + !) + (! + 3w)Y], 


Ym = 2Y[K + X -Y - 


X - 1 
n — 1 


Z,j\zi = 2Z 




n — 1 


iX-l) 


K,m = 2K{1 + K -Y). 


(50) 

(51) 

(52) 

(53) 


The fixed points are given in Table IIIII with their associated solutions and stability. Note that there is a one to one 
correspondence to the solutions found in the previous section. In place of these points, there is now one additional 
fixed point, V, which represents a static universe (since W = 0, in these points H = 0^). Since the point V is 
characterized by Z = 0, the only type of Einstein static solution possible in this context is the one in which the theory 
is effectively coincident with General Relativity (with the difference that in this solution k = 0). The presence of V 
and its features are consistent with the results of In particular, we see that the static solutions can, differently 
from the case of General Relativity, be characterized by an non closed geometry. 

The point 2? is unstable for any values of the parameter n, but it is worth looking at the details of its stability in a 
more careful way. As we have mentioned in the previous section, the stability in terms of the barred variables presents 
some differences with respect to the one in the non barred variables. One of the reasons behind this difference is that 
when X < 0, effectively TL (x —H so that the dimensionless time A4 goes backwards with respect to A4. Remarkably, 
however, (and as expected on the physical point of view) such change does not affect the general structure of the 
phase space in the sense that the stability of the fixed points, and particularly the attractors, is generally not affected 
by the properties of the new time parametrisation. 

The exception is for point 'D in the interval 0 < n < I. Since in this interval this point is a repeller the correction 
for the time parameter would make ti an attractor for X < 0 and a repeller for X > 0. The consequence fo this result 
is that orbits will cross the X < 0 through the point V. Therefore in this specific case the presence of V implies the 
presence of bounces, at least in the neighbourhood of this fixed point. 


^ To be completely precise, H = 0 does not necessarily require X = 0. In fact, X can be finite at the bounce if'H = 0ori<’ = 0- 
such solutions have indeed been found to exist in the Palatini-type (f2.4 —f 0) theories l27ll . However in the first case the IIIHI I implies 
F = const, which means that the Palatini term is linear. We have already seen that this case is degenerate for the dynamical system 
point of view. In the second case {F = 0), since the metrics g and g are not related anymore, the formulation of the theory used in our 
approach looses meaning 













11 


Table 111: The fixed points, their associated solutions and stability for / = analyzed with the barred variables. 
Here A stays for attractor, R for repeller, S for saddle,for attractive focus and F/j for repulsive focus. 


Point 


(X,V,Z,n,K) 


Scale Factor 


Energy Density Stability 


A 


B 


|l,2, ^,0,o| a = aoexp [A(t - to)] /r = 0 Fa 

f 1 , Inn 1 1 = l<n<2 

a = ao(t-to) e- p V 0<w<l 


S otherwise 




2>n{w-\-l)—A 


Q _3„(^+l)+3„+5 , 3n(^+l)-3^-5 ’ a = ao {t - fj. = flQ {t - to) ^ S 

(-3ra(iu+l)+3™+5)^ 


V 


0 , 0 , 0 } 


a = ao 


fj, — 0 


R 0 < n < 1 
S otherwise 


C. The hybrid HS model in the barred variables 


Let us analyze then the case of the hybrid HS model. The system in this case reads 


X^M=X 
Ym = ‘^Y 


K -Y + 1- -{l + iw)X 


--{K{l + 2,w)+^l + w)Z 


2Z- 2ar 


Z 


{K + X -Y + 1) 


k,M = 2k{l+K-Y). 


Z{X-1) - 
2Z - 2aY 


^^[(1 - 3w){K + 1) + (1 + 3w)F]}, 


(54) 


As in the previous section, the new system admits a set of fixed points that correspond to the ones of the formulation 
in terms of the variables m (see Table and,as before, in spite of he differences in the behaviour of the time 
variable, no differences in the stability of the fixed points arises. 

In this case the Einstein static solution is represented by point £ which is characterized by a flat geometry and it 
is always a saddle. This suggests that the phenomenon observed in the previous example, i.e. the fact that bouncing 
orbits will pass through the Einstein point, is not present here, and the study of bouncing phenomenology requires 
numerical tools. 

Finally, there is no global attractor which characterized by de Sitter solution. This confirms that the analysis of 
the behaviour of the action made in the beginning of section IHI Cl cannot be valid in general. 


V. CONCLUSIONS 

In this paper we have applied a dynamical system analysis to investigate the cosmological features of hybrid metric- 
Palatini theories of gravitation. Making full use of the relation between the metrics and g^i,, it is possible to write 
the cosmological equation in a particularly compact way. Using these equations it is relatively straightforward to set 
up a dynamical systems approach for a wide range of hybrid metric-Palatini models. In the present paper, we have 
focused our work on two types of functions /: the first is a generic power of the Palatini curvature TZ] the second has 
a form that is inspired by the Hu-Sawicki model for metric f{R) gravity. 

In the first case, the phase space presents hxed points that are associated to solutions already found, for example in 
[^, with a different approach. The dynamical system analysis, however, clarifies their stability and gives information 
of the global behaviour of the cosmology. For example, it appears clear that the cosmology always presents an 
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Table IV: The fixed points and the solutions of the case with / = a in the barred varibales. The value of the 
parameter A is determined by (1551) and the coordinate of the fixed points. Here A stays for attractor, R for repeller, 

S for saddle, F^i for attractive focus and F/j for repulsive focus. 


Point 

{X,Y,Z,n,K) 

Scale Factor 

Energy Density 

Stability 

A 

[l, 2 , 2 a- 2 ^a(Q-bl), 0 , 0 ] 

a = ao exp ]A(t — to)] 

fi = 0 

Fa 0 < a < 49/32 
A otherwise 

B 

[l, 2 , 2 a-k 2 Vn(a-bl), 0 , 0 ] 

a = ao exp [X{t — to)] 

= 0 

Fa 

C 

fi 8 0 0 -^1 

[ 3 ? 95^5 af, gj 

a = ao (t — to) 

M = /ro(t-to)-"^^+“'^ 

S 

V 

r 1 3 w +7 p 1 nl 

[ 31 U-I -4 ’ 6™-(-8 ’ (3m+4)^ ’ 

a = ao {t — to) 

/i = ^o(t - to)“®^^+™^ 

S 

£ 

[ 0 ,i, 0 , 0 , 0 ] 

a = ao 

fj, = 0 

S 


attractor (yl) associated to an exponential solution. In terms of stability, it is worth remembering that A is not a 
global attractor for the phase space, because of the structure of the dynamical system. 

In the second case, the new formulation of the cosmological equations allows the deduction of some interesting 
properties of the model and in particular provide bounds for the values of the parameters. In terms of the phase 
space, we have now two attractors related to de Sitter solutions. However those fixed points are not global attractors 
for the cosmology. In this respect, the dynamical system analysis reveals the limitations that conclusions drawn of 
the behaviour of the action for small and large values of the Ricci scalar. 

Both the examples we considered present fixed points that correspond to the GR limits of these models i.e. are 
characterised hy Z = 0 (/ = 0). Such points are never stable, so in principle there is no value of the parameters and 
set of initial conditions which might lead to cosmic histories in which the cosmology becomes indistinguishable form 
GR. This is in fact a desired behaviour as we aim at obtaining cosmologies which are close to GR at early times and 
at late time depart from Einstein cosmology. 

The analysis in the alternative ’’barred variables” allows the study of the stability of the Einstein static universe in 
the context of these theories. It turns out that both the models we have analysed contain Einstein static fixed points, 
which are unstable. Their stability however is complicated by the choice of the time variable Ai: in the case V < 0 
this parameter has an opposite sign with respect to Af and the stability analysis must be corrected accordingly. In the 
case of / = this leads to an interesting phenomenon in terms of bounce phenomenology: for a specific interval of 
n orbits close to the Einstein static fixed point cross the V = 0 line only at this fixed points. This result is important 
in terms for the detailed study of bounce dynamics and will be the focus of future studies. 
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